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We performed NLO QCD analysis of the nonsinglet part of the combined SLAC/BCDMS data 
on i*2 with the extraction of a s and high twist contribution. It was shown that the value of a 3 
obtained in the analysis is sensitive to the statistical inference procedures dealing with systematic 
errors on the data. The fit with the complete account of point-to-point correlations of the data gave 
the value of a s (Mz) ~ 0.1180 ± 0.0017(68%C.L.), to be compared with the previously reported 
value of a 3 (M z ) = 0.113 ± 0.003(99%C.L.). This new value of a s is compatible with the LEP 
measurements and the world average. The high twist contribution being strongly anti-correlated 
OO with the value of a s became lower than it was previously reported. 
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I. INTRODUCTION 



O 

Q 

It is well known that the value of the strong coupling constant a s (Mz) measured at LEP is larger than the value 
■ of a s (Mz) obtained from the evolution of results of the analysis of the combined SLAC/BCDMS DIS data on proton 
! and deuterium Q laying at lower Q 2 [Q. This discrepancy caused a lot of discussions (see e.g. j3j) and is often 
attributed to the existence of new fundamental particles, which can change the dependence of a s on Q 2 . Meanwhile, 
the value of a s from J| is strongly correlated with the value of simultaneously fitted high twist (HT) contribution. 
This correlation is inevitable if one does not make a sufficient Q 2 cut of data, otherwise the power corrections can 
essentially, if not completely, imitate the logarithmic scaling violation [Q] . The separation of the power and logarithmic 
behavior is complicated in the case of SLAC/BCDMS data analysis because these data do not practically overlap 
On ' and exhibit significant discrepancies in the vicinity of the overlapping regions. To achieve a satisfactory description 
of the data, one is to invent a method to interpret these discrepancies, which is obviously cannot be done without 
some adoptions. The larger is the correlation of the fitted parameters the more sensitive their values are to the 
perturbations of other inputs to the fit and hence any adoption made in the analysis should be accurately clarified. 
The analysis is not absolutely rigorous in the points concerning the inference of systematic errors. The number of 
independent systematic errors for the combined SLAC /BCDMS data set is about 40 and the authors of g combined 
most of them in a quadrature claiming that this would not distort the results. In the present work we investigated 
the effect of this adoption on the bias of the fitted parameters. 



C3 1 II. THE DATA AND THEIR SYSTEMATIC ERRORS 



We analysed essentially the same data set |3|,|3| as in |2j with the minor differences: 

• we used the data on cross sections separated by the beam energies instead of merged data on F2 . For the SLAC 
data we withdraw the merging errors in this way. The BCDMS data within this approach were reduced to the 
value of R = <jl /&t [[§ common to the SLAC data. 

• we imposed the more stringent cut x > 0.3 to prevent additional uncertainties due to a poorly known gluon 
distribution. This cut leaves the data which can be in good approximation described by the pure nonsinglet 
structure functions, which essentially reduces the number of the fitted parameters. At the same time the value 
of a s in the fit to the combined SLAC/BCDMS data is basically determined by the high-x points and we did 
not loose statistical significance of the analysis as one can see from the final results. The cut x < 0.75 coinciding 
with [^| and rejecting the region where the binding effects in deuterium can be important was also imposed in 
the analysis. The Q range of the data left after the cut is 1 GeV 2 < Q 2 < 230 GeV 2 . 

The number of data points (NDP) and the number of independent systematic errors (NSE) for each experiment 
used in our analysis are presented in Table I. The systematic errors on the BCDMS data are presented by the 
following independent sources: calibration of the measurement of the incident and scattered muon energy, resolution 
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of the spectrometer, detector and trigger inefficiencies, relative normalization of data from internal and external 
targets, general normalization and relative normalization uncertainties between the data set taken at different beam 
energies. (The latest were ascribed to the data at beam energies of 100, 120 and 280 GeV while the data at 200 
GeV were considered as the reference ones.) In the analysis M the systematic errors from the first three sources were 
combined in quadrature into a single error called a "main systematic error" and the data points were shifted by the 
value proportional to this combination while the proportionality coefficient was determined from \ 2 minimization. 
The general normalization was also considered as a free parameter and then the value of normalization uncertainty 
presented in the source paper || was not explicitly accounted for. The rest systematic errors were considered as 
uncorrelated and were combined in quadrature with statistical errors. 

The correlated systematic errors on the SLAC data arose due to: background contamination, spectrometer accep- 
tance uncertainties and radiative corrections uncertainties. In addition, as far the older SLAC data were normalized 
to the data from the E-140 experiment, there are two more systematic errors on them: target dependent and target 
independent relative normalization uncertainties. (The data from E-140 experiment have only one additional absolute 
normalization error). In the analysis 0] all these errors were combined in quadrature with statistical ones. 



III. FITTED FORMULA 



The QCD input leading twist (LT) structure functions of proton and neutron were parametrized at the starting 
value of Qq — 9GeV 2 as follows: 



F£(x,Q ) = A p x a -{l-x) b - — 



F?(x,Q ) = A n x a "(l~x) b " — . 



Here conventional normalization factors N p and N n are 



dxx a ^- x {\ - x) b ?-< 



These distributions were evolved through the region of Q 2 occupied by the data in NLO QCD approximation within 
A/5 factorization scheme with the help of the code used earlier ||. The Q 2 dependence of a s was calculated as 
the numerical solution of the equation 
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a s (Q) a s (M z ) 2tt \M z 



■/?ln 



/3 + l/a s (Q) 



where 



1% 



11 - -n f , (3 



f3 + l/a s (M z ) 
2ttA) 



(3.1) 



51 - ^-n 



and the number of the active fermions n t was changed at the values of Q equal to quark masses keeping the continuity 
of a s . The final formula for structure function used in the fit with account of twist-4 contribution was choose the 



same as m 



F, 



(P.D).HT j-,(P.D).LT 
— ro 



1 , h^\x) 



where F^ P < D '' LT are the leading twist terms with account of the target mass correction ||. The functions hS P ' D \x) 
were parametrized in the model independent way: their values at x = 0.3, 0.4, 0, 5, 0.6, 0.7, 0.8 were fitted, between 
these points the functions were linearly interpolated. As we mentioned before, we used the common value of R (^] for 
all the data including BCDMS ones. 



1 We checked that extra polynomial- type factors do not improve the quality of the fits. 
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IV. RESULTS 



A. BCDMS reanalysis 



At the first stage of our analysis we used the inference procedures analogous to 
through minimization of the functional 



x 2 = E 

KA 



[(/, - X K Ay t )/^ K 



y% 



The parameters were evaluated 



(4.1) 



where K runs through the data subsets obtained by separation of all analysed data on experiments and targets; i - 
through data points within these subsets. The other notations are: yi - the measurements, Cj - the statistical errors, 
combined with some systematic errors as described above, fi - theoretical model prediction depending on the fitted 
parameters, Ay^ - the "main systematic error" on the BCDMS data, \k and are fixed at 0. and 1., correspondingly 
for the SLAC experiments and are the fitted parameters for BCDMS. For the test purposes we fitted formula with 
the parameters £ and A fixed at their values as given in || . The obtained results are presented in column 1 of Table 
II and on Fig. 1. The values of HT coefficients obtained in the analysis [|| are also presented on Fig. 1. As far the 
errors quoted for them in correspond to the change of \ 2 equal to 9., their pictured errors are scaled by the factor 
of 1/3 to provide a meaningful comparison with our figures. One can see that they coincide within the statistical 
fluctuations. 

The next step was to release these parameters (the results are presented in column 2 of Table II). We can note 
that for this fit the BCDMS data are renormalized slightly smaller. As a consequence, the value of a s , which exhibits 
negative correlation with this normalization factor became slightly less than that in || . In this connection note that 
one could suppose the dependence of the normalization factor on the x-cut because the x-shape of the BCDMS data 
does not match the SLAC data very well (in particular, it was pointed in pl[). The errors of the a s value increased 
two times comparing with the first fit. This is in accordance with the above observation, that a s is strongly correlated 
with the normalization factors for the BCDMS data - releasing the latest we allowed more room for the a s variation. 

An alternative possibility to account for the normalization error of the data is to introduce the correlation matrix 

Cij — 5ij(Ti<7j + fifjSj( 



into the minimized functional in the following way: 



K,i,j 



(4.2) 



where sk is the data normalization uncertainty for each target as it is estimated by the experimentalists and Eij is 
the inverse of CV,-; j runs through the data points of each data subset, is the Kronecker symbol and the other 



notations are the same as in (4.1). This approach is natural if one considers a systematic error as a random variable, 
i.e. within the Bayesian approach (see more in ]To| on this scope). The fit within this approach is, in principle, more 
stable comparing with the renormalization approach fl4.1| ) because in (4.2) the normalization parameter variation 
is limited by the scale of s. In our particular case this is not so important as far one can see from Table II, that 
the normalization factors for the BCDMS data are anyway within their norm alization systematic error (3%). This 
anticipation is supported by the results of the fit within the approach (4.2) which are also presented in Table II 



(column 3). Analogously the fitted parameters should not be sensitive to the the stabilization term (£ — l) 2 /s 2 added 



to functional ( LI ) in H as far this term corresponds effectively to the additional measurement of £ with the average 
of 1. and the error of s; the weighted averaging of this measurement with the value of £ from Table II cannot evidently 
change the latest one. 

To proceed with the implementation of Bayesian approach for the treatment of systematic errors, we minimized 
the functional 



X 2 = ^2(fi-Vi)Eij(fj-Vj), 

K,i,j 



(4.3) 



where is the inverse of the correlation matrix 



C > 



SijCTiCTj + fifj(sf ■ sf) 
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and each 4-component vector sf includes the normalization uncertainty as well as the three systematic errors which 
were initially combined into the "main systematic error" of the BCDMS data. The most interesting difference of this 
fit results (presented in Table II, column 4) from the previous fits is the increase of a s . The value of a s is strongly 
anticorrelated with the HT contribution at high x and naturally the last-named decreases correspondingly. The effect 
is of the order of one standard deviation (as could be anticipated because the value of A is of the order of 1. when it 
is released in the fit), with the tendency to decrease the discrepancy with the LEP data. Alongside one can observe 
the decrease of x 2 , which is connected with the fact that in the earlier fits main systematic errors were, as a whole, 
underestimated when combined in quadrature. 

An additional improvement is to account, within this approach, for two more BCDMS systematic errors, which 
were not included in the "main systematics" : The errors due to detector and trigger inefficiencies. The results of this 
fit are presented in column 5 of Table II. Again we can see the enlargement of a s value and the correlated decrease 
of the HT contribution, although not so large as in the case of the re-account of "main systematics" . 

The next step of our analysis was to re-account the errors corresponding to the uncertainty in the relative normal- 
izations of the data subsets for different energies. The results are presented in column 6 of Table II. The value of 
as again increased and the effect is even more pronounced than in the case with the re-account of "main systemat- 
ics" . This is not surprising because as was stated by the BCDMS collaboration itself the uncertainty in the relative 
normalizations have the most effect on the error of a s flijfl . 

Our final exercise with the BCDMS data concerns the correlation of systematic errors on the data from the proton 
and deuterium targets. The authors note that this correlation is large, but do not quantify it. To investigate the 
scale of this correlation effect, wc performed one more fit assuming the total correlation (column 7 of Table II). The 
parameter estimates for real proton/deuterium correlation lie between the values from column 6 and 7, more close to 
7 and we again observe the increase of a s . 

Summarazing, we can conclude that a complete account of point-to-point correlations due to systematic errors on 
the BCDMS data in the combined SLAC/BCDMS analysis cancels the discrepancy with the LEP results. The effect 
of a s increase comparing with the previous analysis |3 arises mainly due to re-account of "main systematics" and the 
errors due to relative normalizations of the data taken at different energies. 



B. SLAC reanalysis 



For the completeness we accounted for the point-to-point correlation of the SLAC data too. At first we proceeded 
with the systematic errors on the E-140 data only. The results of the fit are presented in column I of Table III and 
do not essentially differ from the previous fit. As mentioned above the older SLAC data were renormalized to the 
data from E-140 experiment |Tj|j . Due to the absence of E-140 proton data the renormalization of proton data subsets 
was performed using "bridging" through the E-49B experiment, which introduced additional uncertainties. As far we 
used more of the proton data in the analysis, we preferred to perform the independent renormalization. Then, wc 
removed from the systematic errors on the older SLAC data the relative normalization uncertainties which arose due 
to their renormaliza tion to E-140 and introduced the fitted normalization parameters for each experiment and target 



into the functional (4.3) 



X 2 = E (fifa - Vi) E iAf3^K -Vj), 

K,i,j 

where are fixed at 1. for the BCDMS and E-140 data subsets. The results of this fit are presented in Table 
III, column 2. One can see that our renormalization factors are, as a whole, compatible with 1. within the errors, 
although there is some tendency to shift proton data up comparing with |Q. 

The final step of our analysis was the incorporation of the rest systematic errors into the correlation matrix. 
The results of this fit are presented in column 3, Tabic III. The value of a s due to the last improvement remained 
unchanged, the main effect was a certain increase of % 2 , while the statistical confidence of the fit remains good. This 
is readily understood because if one combines the correlated errors in quadratures, the x 2 is underestimated. In 
the final fit the relative normalization of SLAC data is in the range of few percent up comparing with the BCDMS 
data. In the global fits the SLAC data are often used as the reference ones and the BCDMS data are renormalized 
to them and usually are shifted down by few percent. Our renormalization scheme is in principle compatible with 
the commonly used one, except for the general normalization. This discrepancy cannot be clarified if one uses in the 
analysis only the data on DIS as far it is well known that they cannot define the absolute normalization parameters 
very well, moreover, we applied the cut on x in the analysis. Anyway, it is obvious, that the ambiguity in the general 
absolute normalization cannot affect determination of a slope on Q 2 and, hence, change the value of a s . 
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V. SUMMARY 



The final value of a s (Mz) obtained in our analysis is presented in column 3 of Table III: 

a s (M z ) = 0.1180 ± 0.0017(stat + syst). 

It is compatible with the values obtained in the LEP experiments [|| and in the analysis of CCFR data on F3 JlJ] 
with the extraction of HT contribution |l5| , but is in certain contradiction with the results of (|] : 

a s (M z ) = 0.113 ± 0.003(stat + syst). 

For a meaningful comparison it is worth to remind that in the last result the error corresponds to the change of x 2 = 9, 
i.e. three standard deviations and, consequently, the distance between our result and Q is about 3 standard deviations. 
The statistical confidence of our final fit ( X 2 /NDP = 1179/1183) is perfect, while in @ x 2 /DOF = 599/687 and 
hence the value of \ 2 is by more than two standard deviations lower than its supposed mean. This is yet within 
possible statistical fluctuations, but nevertheless can signal about underestimation of \ 2 due to the combining of 
systematic errors in quadrature. 



In |Q] the value of a s (50 GeV) = 0.180 ± 0.008 with the help of the approximate solution of (3.1) 

2tt ln(21n(Q/A) 



a,(Q) 277 



was transformed into the value of 



/3 ln(Q/A)L f3 f3 HQ /A) 



A^ = 263 ± 42 (stat + syst) MeV. 



Our value of a s (bQGeV) — 0.1935 ± 0.0048 can be analogously transformed into 

A^L = 337 ± 28(stat + syst) MeV. 

The correlation matrix of the fitted parameters is presented in Table IV. One can see from the table that the 
correlation of a s with the HT coefficients is rather large. This supports our initial statement that the separation of 
logarithmic and power effects in a scaling violation is unstable under various assumptions. Other effects, not taken into 
account before (e.g. nuclear effects in deuterium), should be investigated before one can elaborate reliable estimate 
of a s from the analysis of these data. 

As far the HT contribution and the value of a s are strongly anticorrelated, the increase of a s , which we observed 
above, is accompanied by the decrease of HT 0. The total effect on the HT magnitude is about a factor of 3/4, 
comparing with the results ||] . In this connection it is interesting to compare our results with the predictions of the 
infrared renormalon (IRR) model |l2|,[l3| ■ This model is known to reasonably reproduce the shape of HT contribution 
obtained in but the absolute value prediction is about 2.5 times higher than the data (see (l3)). The comparison 
of our results with the IRR model predictions is presented on Fig. 2. The model calculations were made in the 
nonsinglet approximation using the structure functions and the value of A^— obtained in our analysis: 

h(x) = pLT^ x Q) J dzC 2 (z)F 2 LT (x/z,Q) 

C 2 (z) = - 4 . + 2(2 + z + 6z 2 ) - 95(1 - z) - 5'{l - z) 
(1 - z )+ 

2 f3 Q 1 MSi 



2 We omitted the correlation coefficients corresponding to the normalization parameters of the SLAC data due to space 
limitation. The full correlation matrix can be obtained from the author on the request. 
3 This effect was also recently observed in the analysis JIt]], where a B (Mz) was fixed at the value of 0.120 
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where Q 2 = 9 GeV 2 , Cf = 4/3, C = —5/3 and F^ T does not include target mass corrections here. One can see the 
improved agreement of the IRR model predictions with the data at x = 0.5 — 0.7. 

In conclusion, the separation of the logarithmic and power scaling violation effects in the analysis of deep inelastic 
scattering data is unstable due to a high correlation of these effects in the Q 2 region where they are both not small. 
The complete account of point-to-point correlations of the data lead to the shift in the value of a a by about 3 standard 
deviations comparing with the simplified statistical inference procedure. The HT contribution, which is strongly anti- 
correlated with a s , decreases within this approach and becomes more compatible with the prediction of IRR model 
at moderate x. Further investigation of a possible perturbation in the analysis of DIS data is needed before a reliable 
value of a s can be determined. 
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FIG. 1. The high-twist contributions obtained in the fit with the functional (1) (full circles and lines). The results of the 
analysis [2] are presented for comparison (open circles). 



FIG. 2. The high- twist contributions obtained in our final fit (full circles) and the results of the analysis [2] (open circles). 
The full curves represent the calculations on the IRR model. 
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TABLE I. The number 
sets. 


of data points (NDP) and the number of indepenc 


lent systematic errors (NSE) for the 


i analysed data 


Experiment 


NDP (proton) 




NDP (deuterium) 


NSb 


BCDMS 






223 




162 




9 


E-49A 






47 




47 




5 


E-49B 






109 




102 




5 


E-61 






6 




6 




5 


E-87 






90 




90 




5 


E-89A 






66 




59 




5 


E-89B 






70 




59 




5 


E-139 










16 




5 


E-140 










31 




4 


TOTAL 






611 




572 
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TABLE II. The results of the fits with the various approaches to the treatment of the BCDMS systemat 


ic errors. The 


parameters £ and A describe the renormalization and shift of the BCDMS data, /13, 4, 5,6, 7, 8 


are the fitted values of the HT 


contribution at x = 0.3, 0.4, 0.5, 0.6, 0.7, Oi 


>. For the description of the columns see the text. 








1 


2 


3 


4 


5 


6 


7 


Ap 


0.612 ±0.028 


0.579 ± 0.028 


0.581 ± 0.028 


0.523 ± 0.024 


0.531 ±0.024 


0.531 ± 0.024 


0.519 ±0.022 


a p 


0.642 ± 0.028 


0.689 ± 0.032 


0.685 ± 0.032 


0.748 ± 0.033 


0.736 ± 0.032 


0.734 ± 0.032 


0.748 ± 0.030 


b v 


3.588 ± 0.029 


3.675 ± 0.038 


3.670 ± 0.038 


3.702 ± 0.038 


3.686 ± 0.037 


3.670 ± 0.037 


3.667 ± 0.035 




4.0 ±3.2 


4.0 ±3.5 


3.7 ±3.0 


4.7 ±4.8 


3.4 ±2.6 


4.2 ± 3.8 


4.7 ±4.4 


a„ 


0.14 ±0.11 


0.14 ±0.12 


0.15 ±0.12 


0.12 ± 0.12 


0.16 ±0.12 


0.13 ±0.12 


0.12 ±0.11 


b n 


3.52 ±0.12 


3.52 ±0.14 


3.54 ±0.14 


3.48 ±0.14 


3.52 ±0.14 


3.48 ±0.14 


3.51 ±0.12 


a s {M z ) 0.1141 ±0.0007 


0.1089 ±0.0016 


0.1093 ±0.0016 


0.1119 ±0.0015 


0.1140 ±0.0017 


0.1173 ± 0.0018 0.1188 ±0.0018 


Ap 


1.4 


0.95 ±0.13 


0.97 ±0.13 










A D 


1.2 


0.89 ±0.15 


0.90 ±0.15 












0.99 


1.0138 ±0.0059 














1.004 


1.0261 ±0.0063 












hi 


-0.154 ±0.016 


-0.136 ±0.017 


-0.138 ±0.016 


-0.114 ±0.017 


-0.125 ±0.018 


-0.136 ±0.018 - 


-0.136 ±0.017 


hi 


-0.009 ±0.019 


0.030 ± 0.022 


0.026 ± 0.022 


0.015 ±0.022 


-0.010 ±0.024 


-0.047 ± 0.026 - 


-0.068 ± 0.026 


hi 


0.175 ±0.029 


0.257 ±0.038 


0.250 ± 0.037 


0.191 ±0.038 


0.149 ±0.041 


0.077 ± 0.045 


0.029 ± 0.046 


hi 


0.623 ±0.054 


0.803 ± 0.072 


0.788 ± 0.070 


0.643 ± 0.071 


0.572 ± 0.077 


0.440 ± 0.083 


0.338 ± 0.084 


hi 


1.106 ±0.089 


1.49 ± 0.13 


1.46 ±0.13 


1.23 ± 0.13 


1.11 ±0.13 


0.90 ± 0.14 


0.73 ±0.14 


hi 


1.83 ±0.25 


2.56 ±0.31 


2.51 ±0.31 


2.20 ± 0.30 


1.99 ±0.30 


1.66 ±0.31 


1.41 ±0.30 


hi 


-0.130 ±0.018 


-0.102 ±0.019 


-0.103 ±0.019 


-0.094 ±0.019 


-0.102 ±0.020 


-0.123 ±0.021 - 


-0.129 ±0.021 


hi 


0.048 ±0.017 


0.104 ±0.022 


0.099 ± 0.022 


0.081 ± 0.022 


0.054 ± 0.025 


0.010 ± 0.028 - 


-0.005 ±0.029 


hi 


0.266 ± 0.027 


0.367 ± 0.038 


0.358 ± 0.037 


0.299 ± 0.038 


0.248 ± 0.042 


0.172 ± 0.047 


0.146 ±0.049 


hi 


0.657 ± 0.050 


0.844 ± 0.069 


0.829 ± 0.068 


0.696 ± 0.068 


0.611 ±0.075 


0.480 ± 0.082 


0.445 ± 0.085 


hi 


1.050 ±0.075 


1.38 ±0.11 


1.36 ±0.11 


1.15 ±0.11 


1.03 ±0.12 


0.82 ±0.13 


0.77 ±0.13 


hi 


2.28 ±0.25 


2.96 ±0.31 


2.92 ±0.31 


2.52 ± 0.30 


2.34 ± 0.30 


1.98 ±0.31 


1.94 ±0.31 


x 2 


1090.5 


1067.5 


1068.3 


963.7 


964.3 


973.3 


971.5 
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TABLE III. The results of the fits with various approaches to the treatment of the SLAC systematic 


errors. The pa- 


rameters f describe the renormalization of the SLAC data, 


^3,4,5 f 6,7 f 8 are the fitted values of the HT 


contribution at 


x = 0.3,0.4,0.5,0.6,0.7,0.8 


. For the description of the columns 


see the text. 






-i 
1 


O 

z 


3 


Ap 


a co7 _i_ n noo 
U.oz ( zb U.Uzz 


n k /i c _i_ n hop: 
U.D4D zh U.Uzo 


A CT 1 G 1 A AOO 


a p 


n 700 _|_ n AQA 
U. 1 oo zh U.UoU 


U. fzo zh U.UoU 
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TABLE IV. The correlation matrix for the parameters from the final fit. 
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